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We review the connection between QCD and potential models of quarkonia in 
the framework of effective field theories, with an emphasis on non-perturbativc 
methods such as lattice simulations. The static and heavy quark potentials are 
introduced and the applicability of the non-relativistic and adiabatic approxima- 
tions are discussed. The roles of so-called hybrid potentials are explored and we 
comment on the QCD analogue of the QED Lamb shift. 



1 Potential models 

Long before direct non-perturbative QCD evaluations became feasible poten- 
tial models, some QCD inspired and some not, turned out to be extremely 
useful for an understanding of the masses and decay rates of bottomonium 
(T), charmonium (J/ip) and B c states. The underlying idea is that when the 
quark mass m is bigger than all other bound state energy scales, time scales 
AtQ that characterise the relative heavy quark movement become much big- 
ger than time scales Ai s , associated with the gluon (and sea quark) dynamics, 
too. In this situation the leading Born-Oppenheimer (or adiabatic) approxi- 
mation should be applicable and feedback effects of the moving heavy quarks 
onto the accompanying gluons and sea quarks might be neglected. If in ad- 
dition the typical relative quark velocity is much smaller than the speed of 
light, v = p/m <C 1, the bound state should be governed by a Schrodinger 
Hamiltonian, 

P 2 

H = 2m+ ^— + V phcn (r), Hi/> n = E n ^ n , (1) 

where ijlr = ra/2 denotes the reduced mass and the potential VphenM param- 
eterises the interaction that is mediated by gluons and sea quarks. 

While in lattice QCD calculations spectrum and wave functions can in 
principle be predicted for a given set of quark masses and scale parameter, 
in the potential model approach one is initially confronted with the inverse 
problem of guessing or deriving the potential tphen from the knowledge of 
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some energy levels and decay rates, before any predictions can be made. The 
experimental spin averaged quarkonia levels seem to constrain the potential (up 
to an additive constant) reasonably well for distances 0.2 fm < r < 1 fm to a 
form that can be represented in terms of a funnel (or Cornell) parametrisation,B 

V phcn (r) =ar- ~, (2) 

with e « 0.5 and a string tension a « 425 2 MeV 2 . The advantage of potential 
models over other approaches such as Lattice Gauge Theory or Spectral Sum 
Rules is the relative ease with which properties of radial excitations can be 
derived and the simplicity and transparency of the representation of the bound 
state problem in terms of a nonrelativistic quantum mechanical Hamiltonian. 

If a potential model has anything to do with QCD, the heavy quark po- 
tential VphcnM should be calculable on the lattice and such a determination 
should increase its predictive power. One object that is easily accessible in 
lattice studies is the so-called static potential, V${r). In this article we shall 
address the following questions: does the "heavy quark potential" Vphen ex- 
ist? What is the relation between V^,hcn and Vb? Can relativistic corrections 
systematically be incorporated into Eq. (|])? 

2 QCD static potentials 

Ironically, forgetting everything one knows about quarks and their properties 
like spin and mass seems to be a good strategy for an understanding of heavy 
quark bound state problems. After successfully mastering this first mental 
barrier one has to place an external colour charge Q and anti-charge Q* at a 
distance r into the QCD vacuum. The static potential Vo(r) is then the energy 
difference of this configuration with respect to the vacuum energy. Since we 
are dealing with a quantum field theory, an additive mass renormalisation is 
required, due to the self energy of the point like sources, 

V , P hys(r) = Vo(/i;r) - V se \f((J,), Keif (n) oc /i/ In n (ju — > oo) : (3) 

the static potential Vo(M! r ) diverges as the cut-off on gluon momenta /x is sent 
to infinity while, once a subtraction scheme is specified, Vo jP hys is well defined. 

Vo(n;r) can be computed from Wilson loops, for instance on a lattice 
with spacing a ~ %/ fx. From perturbative QCD one would expect the generic 
Coulomb short distance form, Vo(r) = —e/r, where the effective e = (4/3)a/j(r) 
will logarithmically depend on the distance r while — neglecting sea quarks — 
one might expect a linear rise at large distances, in agreement with the funnel 
parametrisation, Eq. (0). 
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Figure 1: Comparison between the quenched lattice QCD potential (full symbols) and per- 
turbation theory. Uu 

In QCD the potential is the energy of a bound state of gluons and static 
sources. This contrasts the situation in QED where confinement is lacking and 
photons carrying any surplus energy will just be radiated away. In QCD a 
discrete-spectrum of gluonic excitations exists, the so-called hybrid potentials, 
see e.g. Bel. This leads us to conclude that in addition to the distance r another 
scale, AV(r), the gluonic binding energy, has to be considered. In a quenched 
setting AV(r = 0) will be a glueball mass, the lightest of which is about 
1.6 GeV, while in the case of QCD with sea quarks, radiation of a pair of 
pions becomes possible. While the coupling constant a s (r _1 ) — * as r — ► 
0, a g [AV(r)] remains large at all distances, from which one might expect 
perturbation theory to break down when naively applied to the Wilson loop. 

In Fig. [l] we compare the quenchetLQCD potential, calculated on the lattice 
with continuum perturbation theoryCrfl in units r y ss 0.5 fm at short distances, 
r < 0.23 fm. Note that the A-parameter has been determined from an un- 
related lattice observable, such that the only freedom in this representation 
is an additive constant, due to different values of V se \f in different regularisa- 
tion schemes. Indeed, the perturbative series does not seem to convincingly 
approach the non-perturbative result and the coefficients of the perturbative 
expansion are large as indicated by the big difference between one- (order a 2 s ) 
and two- loop (order a 3 s ) potentials. This means that different resummation 
prescriptions will in general yield different two-loop results, some being closer 
to the non-perturbatively determined potential, some disagreeing even more. 
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Figure 2: Static QCD potential and first hybrid excitation, without (quenched) and with 
(k = 0.1575) sea quarks. Q 

Gluonic excitations can be classified with respect to representations of 
the cylindrical symmetry group Z?oo/i- The lowest hybrid excitation, in which 
the glue carries one unit of angular momentum about the intersource axis, is 
labelled by H u . A recent lattice result □ of both the static potential and 
the hybrid potential U u , normalised such that Vb(r ) = is depicted in Fig. |[ 
both without sea quarks (quenched) and with 2 flavours of sea quarks, slightly 
lighter than the strange quark (k — 0.1575). 

When including sea quarks the potential will eventually flatten around i 
r ~ 2.3 ro (string breaking) while in the quenched case the linear rise will 
continue. Within the region r > 0.2 fm the lattice potentials can be fitted to 
the parametrisation Eq. (g), with the result e « 0.3 neglecting sea quarks and 
e w 0.36 with two flavours of light sea quarks. 

3 The infra red problem of perturbation theory 

In QED without fermions, the perturbative expansion of a Wilson loop can 
schematically be written as, 
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multi-photon exchanges exponentiate! The solid lines denote the static sources 
and an integration of the vertices over all positions along the contour of the 
Wilson loop is understood to be included within the notation. The static 
potential in the limit of large Euclidean times r is related to the Wilson loop, 

(W(r, t)) oc exp[— Vo(r)r] (r — > oo), (5) 

which means that only contributions whose logarithms are proportional to t 
have to be considered; interactions with the spatial closures of the loop can be 
ignored. The single photon exchange within the exponent on the right hand 
side of Eq. (^) can easily be calculated: a factor r is obtained by transla- 
tional invariance and the remaining time integral yields a <5-function for the 
4-component of the photon momentum. After a 3-dimensional Fourier trans- 
form one finds the well known result, Vo = —af s /r to all orders of perturbation 
theory. 

The situation becomes somewhat more involved if fermions are included, 
due to the occurrence of diagrams such as, 



D 



(6) 



However, irreducible diagrams Di still exponentiate and the relation, Vo(r) = 
— Hm T ^ 00 '^2 i dDi(r,T)/dT, holds. In a non-Abelian theory, starting from a 
two gluon exchange, the colour pre-factor of a given diagram depends on the 
ordering of the vertices along the contour of the Wilson loop. Hence, contri- 
butions are encountered that do not exponentiate anymore. □ 

Still, at least up to order af, the static potential Vq agrees with the result 
one obtains ignoring gluons that couple to the spatial closures of the Wilson 
loop, the so-called singlet potential V s . Within perturbation theory it is also 
possible to define a so-called octet potential, V , by introducing colour gener- 
ators at zero and infinite times into the colour traces. In doing so, one finds 
V = —V s /8, at least to order a 2 . It is possible to calculate hybrid potentials 
in perturbation theory too, starting from linear combinations of Wilson loops 
where the spatial connections are deformed or colour fields are inserted. While 
for instance the hybrid potential Vn u vanishes to order a s , the next excited 
state, Vr- agrees with V to this order. 

Up to order non-Abelian contributions do not cause any fundamental 
problem, however, it has early been noticedQ that diagrams, 
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diverge faster than r as r — > oo. This means that a nai've perturbative ex- 
pansion of the static potential encounters an infra red problem: between the 
times t\ and ti the quark and anti-quark are in a colour octet state and the 
first argument of the time integral should be proportional to exp(— AVt) with 
AV ~ V — V s , which is small for times t = t% — t\ ^> AV^ 1 . However, 
in perturbation theory the real space propagator for large times only decays 
like l/(r 2 + t 2 ) and this (wrong) behaviour finally causes the infra red diver- 
gence. It has therefore been suggestedJp regulate the encountered divergence 
by resumming all possible interactions!!! between t\ and ti . 

Since the infra red problem is related to the emission of an ultra-soft gluon, 
it has also been equated to the Lamb shift in the literature. □ One should, how- 
ever, keep in mind that in QED the Lamb shift affects atomic bound states and 
not the Coulomb potential itself (which in QCD happens to be a bound state 
too) . We shall see below that another "Lamb shift" , analogous to that of QED 
exists within quarkonia systems. One should also remember that while pertur- 
bation theory generically encounters problems when applied to bound states, 
the static potential is perfectly well defined in a non-perturbative context. 
Despite all problems with the static potential itself, Potential Nonrelativistic 
QCD0 (pNRQCD) might allow for a consistent and systematic perturbative 
treatment of heavy quark bound states. Within the pNRQCD Lagrangian 
a heavy quark singlet potential \/p nr Q cd appears that can be related to the 
static singlet potential V s . It has been demonstrated E2l that the divergence en- 
countered in V s at order is cancelled by a counter term from the matching 
of pNRQCD to QCD: the resulting y" s pNR Q CD does not incorporate any singlet- 
octet-singlet transitions which are dealt with explicitely in the pNRQCD mul- 
tipolc expansion. ._. 

Finally, I mention the remarkable fact that up to order EJ potentials 
between charges within different representations D of the gauge group scale 
in proportion to the respective Casimir factors, Co = (1/ -D)tr T®T® . This 
Casimir scaling seems to be rather accurately satisfied in the non-perturbative 
regime too.E 2 ] 

4 The "heavy quark potential" from QCD 

The natural starting point for deriving a non-relativisticJ^amiltonian for quark- 
onia bound states is an effective field theory, NRQCD:I12 

T D 2 



£nrqcd — — ^ 



iP + xH---]x + £ s iuc (8) 

2m Q J 

ip and x are the quark and antiquark Pauli spinors. We eliminate the matching 
constant C2 by defining m = mg/c2 = tuq + Sm, where vtlq is the bare quark 
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Table 1: Hard, soft and ultra-soft scale estimates. 





J/i/> 


T 


tt 


e + e 


to 


1.4 GeV 


4.7 GeV 


175 GeV 


511 keV 


mv 


0.7 GeV 


1.3 GeV 


45 GeV 


3.7 keV 


mv 2 


0.4 GeV 


0.4 GeV 


12 GeV 


0.027 keV 


V 


0.5 


0.29 


0.26 


0.007 



mass appearing in the original Dirac Lagrangian in Euclidean space, C — 
q{Ip+m,Q)q. The underlying idea of an effective field theory description is that 
if the quark mass is much bigger than exchange momenta and binding energy, 
QCD processes of scales q > fi > mv can be integrated out into local matching 
constants Cj. 

For systems containing only one heavy quark the leading order effective 
Lagrangian will not incorporate the kinetic term and the usual heavy quark 
effective theory (HQET) power counting in powers of the inverse quark mass 
applies. However, heavyonia bound states cannot be addressed within a purely 
static theory, basically because, without a kinetic term, quarks that travel 
at different velocities will never interact with each other: hence, the leading 
order Lagrangian is that of Eq. . Since the static theory does not apply to 
quarkonia systems it is also not a priori clear whether the phenomenological 
potential V^henrAf Eq. (Q) is at all related to the static potential Vq in some 
approximation. U To keep the discussion brief I will use the conventional power 
counting rulesEj in the relative heavy quark velocity v throughout this article. 
Several well founded objections to this scheme exist and ot/hpr /? n ]JT[ting schemes 
have been suggested to be more suitable or better defined. Eju tSElO According 
to the "standard" conventions the leading order NRQCD Lagrangian is counted 
as order v 2 . 

In Table [l] estimates of various scales for the charmonium, bottomonium 
and (unstable and, therefore, hypothetical) toponium ground states from a 
potential calculation u are displayed. For comparison, we include the corre- 
sponding estimates for positronium. While binding energies and level split- 
tings are of size mv 2 (ultra-soft), the momenta exchanged between the quarks 
are of order mv (soft). NRQCD is the effective theory for physics at scales 
q < jj,, mv < /i < to, while physics of scales q > /x, has been integrated out 
into matching coefficients Cj(m//x, a s ) that are in principle calculable from the 
QCD Lagrangian. It is possible to integrate out the soft scale mv in aJiirther 
step with the result of the effective theories, pNRQCDB or vNRQCD.H 

In NRQCD, being a nonrelativistic theory, the time evolution of a quark 
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in a background gauge field, {A ll {x)} 1 is controlled by a Hamiltonian, 

D 2 

Hj, — m — 5m + igA^ h • • • ■ 

2m 



(9) 



The quark propagator, K — ifnp obeys the evolution equation, —d^K = H^K, 
which for the initial condition, K(x, y)\ Xt=Vi = £ 3 (x — y), is solved by, 

K(x,y)oc DzDpexpl dt [pz - jEfy(z,p)] \ , (10) 

where the path integral is over all paths connecting y with x. 

For illustrative purposes we sketch the derivation of a quantum mechanical 
Hamiltonian for quarkonium systems to leading order in v: we combine two 
propagators for quark and antiquark into a four-point function, 

Xi=(0,t) >fe= (r,x) 



Ut 



G(r, 



Uo 



-Retr 



{u^K(y 1 ,x 1 )U T K(x 2 ,y 2 )} (11) 



y = o 



V = (r,0) 



U T denotes a gauge transporter at time r. To higher orders in the velocity 
expansion G will not only depend on the distance but also on spin, angular 
momentum and momentum. It can be shown Ea that to leading order only 
the shortest paths between j/j and Xi contribute to the path integral Eq. (^0|), 
such that after factorising the gauge field dependent and independent parts 
one obtains, 



(G(r, t)) cx exp 



r T _ 

7 dt H 

Jo j=1 



Zj — 5m j 



2m i 



(W(r,r)), (12) 



where W(r, t) is the usual Wilson loop and the expectation value is over gauge 
field configurations. Now the Hamiltonian H that governs the evolution of the 
two-particle state, H(G) = —d4(G) can easily be read off, 



H = 2[m - 5m(jj)] + ^ + ^o(m; r), 

2fJ-R 



(13) 



where we have set m = mi — m 2 — 2//r. Indeed, the leading order Hamil- 
tonian is that of Eq. (|j) with V p hcn(r) — Vb(/i;r) — 25m(n). The quantum 
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Table 2: Combinations of spins and angular momenta that can couple to J c = 1 
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mechanical Hamiltonian should not depend on the scale /i. This implies the 
relation, 2Sm(ljf) — 25m(fi) + V sc \{(n') — V sc u(n), between the mass shift Sm(fi) 
between "kinetic" mass m and QCD quark mass rriQ (/i) and the static self en- 
ergy V scif (fj,), defined in Eq. (§). 



5 The "Lamb shift" and relativistic corrections 

The assumption, At g > Atg, underlying the adiabatic approximation, is 
clearly violated whenever ultra-soft gluons with momenta smaller than the 
interaction energy are exchanged; at the end of an interaction time of order 
AtQ such a gluon will still be present and cannot be integrated out into a 
potential. In the positronium case such effects result in the Lamb shift which 
can be accounted for by considering transitions between higher Fock states, 
|e + e~7), |e + e _ 77), . . . and |e + e~). Any formalism that does not incorporate 
such transitions can only be approximate. Where then is the loop hole of 
the derivation of Sec. || above? The definition of the 4-point Green function 
Eq. (|ll|) is ambiguous since the Pooh states of the gluonic connections Uq and 
U T are unspecified. In general, G will carry two more indices that run over 
gluonic excitations, a,b = S+, IL, • • •, in addition to time, distance, spin, 



momentum and angular momentum of the quarks. The derivation of Eq. (13) 
can then be generalised, resulting in the Schrodinger equationtl H^ip^ = Ei/j®, 
where the J PC quantum numbers of quarkonia eigenstates ip n are products of 
quark and gluon quantum numbers. Physical states will be linear combinations 
of standard quark model states (Sj) and various hybrid excitations. 

The off-diagonal elements P s +& mediate transitions between standard 
states and states with excited glue and are related to gluon emission by the va- 
lence quarks. Each such element is accompanied by a pre-factor gv, g = y/Anos. 
We may therefore treat mixing effects as perturbations. Let us start from the 
unperturbed, diagonal Hamiltonian, 

H aa =2(m- 5m) + V a (r), H aa ^ a = P°'°V°'°. (14) 

The energy shift of the state n of the lowest lying channel (a = £+) is, 
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AF a - \(^n a \Hgb\ 1pm b )\ _ ~ I 2 yO,a\ 

AL n - =o^ =q^ - V [UsV A n ) , (15) 

rn,b^a £ " m £/ " 

where only states with equal J PC content yield a non- vanishing matrix element 

n 

in the nominator and E n = E n ' 9 +AE n 9 +0(v 6 ). In Tableg we have listed all 
possible combinations of gluonic excitation A^ 1 ' , total angular momentum 
K and spin S of a vector particle (J PC = 1 — — ). It turns out that to order 
v only mixing-.of the standard <S* wave state with the S~ P wave state has to 
be considered. U In this formalism spin exotic states like 1 that do not have 
any E+ content are automatically accounted for. 

Having understood these mixing effects, we can now say that "valence 
gluons" that accompany the quarks in the form of hybrid excitations of the 
flux tube and "sea gluons" , whose average effect can be integrated out into an 
interaction potential, can be distinguished from each other. To lowest order of 
the relativistic expansion pure quark model quarkonia and pure quark-gluon 
hybrids exist, which then undergo mixing with each other as higher orders in v 
are incorporated. Moreover, at leading order the Hamiltonian is of Schrodinger 
type and Vphcn = Vq. As soon as one allows for higher orders in the velocity 
expansion, an instantaneous interaction potential does not exist anymore but 
the corrections can systematically be accounted for. 

With respect to QED the Lamb shift is naively enhanced by a relative 
factor, a s v\ / '(ai s v^ + ) w 6 x 10 4 . However, below a certain glueball (or, 
when sea quarks are included, meson) radiation threshold, the spectrum of 



excitations is (unlike in the QCD case) discrete and the denominator of Eq. (15) 
will become large. It is clear that as long as gaps between hybrid potentials arc 
bigger than quarkonia level splittings with respect to the IS* state, transitions 
will be energetically penalised: the wave function will only see the lowest 
lying potential and the adiabatic approximation is reliable while one might 
expect levels such as T(3S') and T(4S I ) to be affected by the presence of hybrid 
excitations. 

We notice formal similarities between the discussion above and the pN- 
RQCD Lagrangian that includes transition elements between two different 
states, singlet and octet. However, to draw the analogy, singlet ~ E^" and 
octet ~ hybrid, is not straight forward. Moreover, in the potential frame- 
work transitions between different -Dqo/i states are only possible via interac- 
tions between the valence quarks and the glue while in pNRQCD an ultra-soft 
gluon that causes a singlet-octet transition can also be radiated from another 
gluon. Effects of this latter kind are already included automatically in the 
non-perturbative matrix elements H a b, defined through expectation values of 
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Wilson loop like operators. 

The velocity will approach zero only as a logarithmic function of the quark 
mass: the numbers of Table [j] reveal that v 2 is reduced from 0.085 to just 0.070 
when replacing bottom quarks by (almost forty times heavier) top quarks. 
Therefore, although v — > as m — > oo, at any reasonable quark mass value one 
has to address Lamb shift effects as well as standard relativistic corrections 
which are also mandatory for an understanding of fine structure splittings. 

The complete Hamiltonian to order i> 4 is knownieiiffl where var- 
ious corrections are given in terms of expectation values of W-ilson loop like 
operators, most of which have been calculated on the lattice. c3l3 In addition 
to Pauli, Thomas and Darwin like terms, momentum dependent corrections 
appear due to insertions of the operators D4 and D 2 /(2m) as well as 1/m and 
1/to 2 corrections that are spin- and momentum-independent. El We remark 
that the momentum dependent corrections have nothing to do with thq-X.amh 
shift. The derivation can either be performed alofts-thc lines of Sec. ||lJ'E3ll3 
or in terms of quantum mechanical perturbations . H'EoO The predictive power 
is at present in particular limited by an insufficient knowledge of the matching 
constants between QCD and NRQCD. For instance we estimate E3 a related 
uncertainty of 25 % on T fine structure splittings while relativistic order v 6 
corrections are unlikely to change order v 4 QCD potential predictions by more 
than 10 %. 



6 Summary 

We have demonstrated that in a non-relativistic situation, it is possible to 
factorise gluonic effects from the slower dynamics of the quarks and to derive 
"potential models" from QCD. Effective field theory methods turn out to be 
essential for this step. The resulting Hamiltonian representation of the bound 
state problem in terms of functions of canonical variables offers a very intuitive 
and transparent representation of quarkonium physics. It highlights parallels 
as well as differences to well understood atomic physics. 

The adiabatic approximation is violated when ultra-soft gluons are radi- 
ated, i.e. when the nature of the bound state changes during the interaction 
time. Such effects can systematically be incorporated into the potential formu- 
lation by enlarging the basis of states onto which the Hamiltonian acts and a 
lattice study of the relevant transition matrix elements H a b is on its way. The 
validity of the adiabatic approximation is tied to that of the non-relativistic 
approximation in so far as off-diagonal entries of H a b are suppressed by pow- 
ers of the velocity, v. Hybrid mesons become a well defined concept in the 
potential approach and translation into the variables used for instance in flux 
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tube models is straight forward. In view of phcnomcnological applications, a 
precise determination of the matching coefficients between QCD and (lattice) 
NRQCD is required. 
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